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Large scale motion of cytoplasm called cytoplasmic streaming occurs in some large eukaryotic 
cells to stir the cells constituents. In Drosophila oocytes, microtubules have been observed to 
undergo undulating motion, curving to form travelling waves during cytoplasmic streaming. Here 
we show that this wave- like motion can be understood physically as due to the hydrodynamic drag of 
streaming impellers attached to kinesin motors moving toward the plus-ends of microtubules whose 
minus ends are anchored to the cell cortex. The tangential forces applied to such microtubules 
by kinesin give rise to bending and leads to chiral symmetry breaking causing the microtubules to 
propagate long travelling waves. The waves are reminiscent of those seen in flagellar motion but of a 
much longer time scale and by a different physical mechanism. We show how kinesin movement can 
produce a bulk flow of cytoplasm surrounding a microtubule with the range of flow greatly enhanced 
by the effect of hydrodynamic coupling between impellers. That is, a relatively small number of 
motors can move a large amount of fluid. The chaotic nature of the fluid motion of cytoplasm caused 
by kinesin movement along constantly changing microtubule trajectories is important as it greatly 
enhances the efficiency of mixing. Existing data on in vitro microtubule gliding assays also show 
this chiral instability in two dimensions and an analysis of this gives quantitative estimates for the 
forces exerted by motors and the drag coefficient. 



I. INTRODUCTION 

Microtubules are flexible hollow polymers of tubulin 
subunits that serve many critical functions in eukaryotic 
cells. They are utilized in structural contexts, because of 
their relatively stiff, yet flexible, mechanical properties. 
They also act as directional highways through the vis- 
cous cytoplasm. Molecular motor proteins carry cellular 
constituents along the microtubules with kinesin moving 
toward their fast growing "plus-ends" and dynein moving 
toward their slow growing "minus-ends" . Many studies 
have focused on motor driven transport processes that 
generate asymmetric distributions of specific cytoplasmic 
constituents; asymmetries that are essential for complex 
cellular functions. In this letter we will analyze a surpris- 
ing role for microtubules and kinesin in a mass transport 
process called cytoplasmic streaming that has evolved to 
accomplish just the opposite; efficient homogeneous mix- 
ing of the contents of a cell [l| and see there Supplemental 
Movie 13 0. 

Vigorous streaming is initiated during the final stages 
of development of Drosophila oocytes to disperse asym- 
metrically distributed mRNA particles, protein com- 
plexes, and membranous organelles. The mixing process 
is important for subsequent embryonic development, and 
cannot be accomplished by diffusion alone. For example 
a 1/im yolk-filled vesicle in the cytoplasm, assuming a 
viscosity 8 times that of water [3| would take approxi- 
mately a week to diffuse the 500/xm length of an oocyte. 
This is far too long to satisfy the need for the mixing 
of yolk-filled oocyte cytoplasm with the mass of yolkless 
nurse cell cytoplasm that floods the anterior of the oocyte 
near the end of its development. 



The problem of mixing in small systems, such as in 
microfluidics chambers, has been the subject of much in- 
vestigation 3. The way that a fluid at low Reynolds 
number is stirred has a profound effect on the efficiency of 
homogenization. For example, the steady state fiow fields 
generated by a single stir bar inside a closed chamber are 
far less efficient than the more chaotic flows generated 
when several stir bars are used 0, . Rigorous analyses 
in two dimensions show that mixing is most efficient when 
topological chaos is created by three or more stirrers, 
as can be shown by application of the Thurston-Nielsen 
classification theorem With this in mind, it is in- 

teresting to note that during fast cytoplasmic streaming 
in Drosophila oocytes, microtubules appear to be locally 
aligned along dynamically changing curved paths that 
produce travelling waves [l[ . The streaming cytoplasmic 
fiuid moves along those paths, in patterns reminiscent of 
fiowing water and seaweed. Yolk particles in the cyto- 
plasm have a speed of roughly 0.25^m/s. The particles 
within a region stream for the most part in one direc- 
tion, but with a non-negligible deviation in that direc- 
tion over time. The fluctuating directions, which parallel 
the curved paths of the microtubules in the same region, 
serve to stir cytoplasm in a chaotic manner that, as the 
preceding paragraph suggests, is important for efficient 
mixing. 

At first sight, it might appear that the time-dependent 
wave-like motion is due to turbulence of the surround- 
ing fiuid. However at such minuscule Reynolds numbers, 
inertial effects are negligible and turbulence is impos- 
sible Therefore one is left with a mystery of the 
relationship of fluid and filament and how such chaotic 
patterns could come about. The plus-end directed mo- 
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tor kinesin-1 plays a crucial role in this, as shown by 
genetic mutations in its force producing subunit (Khc) 
that prevent streaming and mixing The speed of 
unloaded kinesin along the microtubule has been mea- 
sured to be in the 0.5 — 1/im/s range [ll|, [T^]- This 
is higher than the fluid speeds measured during fast cy- 
toplasmic streaming [l| and is also consistent with the 
role of kinesin in powering the mixing. Inhibition of the 
opposing minus-end directed motor protein, dynein, has 
a complementary effect, actually stimulating fast cyto- 
plasmic streaming Microscopy studies suggest that 
microtubules participating in this motion have their mi- 
nus ends attached to the cortex and their plus ends away 
from the cortex in the interior of the oocyte [H, [l^l . 

The explanation that we analyze for the streaming 
phenomena is very simple: the mass motions of cyto- 
plasm and the microtubule undulations are complemen- 
tary physical consequences of kinesin moving cargo to- 
ward the plus ends of microtubules whose minus ends 
are in contact with the cortex. The cargoes serve as im- 
pellers that both drive the fluid motion away from minus- 
ends and generate tangential forces that move plus-ends 
toward minus-ends causing bends in the microtubules. 
This has been suggested previously to explain cytoplas- 
mic streaming, but without a physical model of the mech- 
anism [l| . The analysis below shows that long range hy- 
drodynamic forces couple individual impellers, resulting 
in an effective mechanism for drag-induced bulk move- 
ment of cytoplasm. We then show that an instability in 
the dynamics leads to chiral symmetry breaking giving 
rise to wave-like motion of microtubules, and show that 
the time and length scales predicted are in good agree- 
ment with the previous experimental results. 
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FIG. 1. A train of spherically shaped impellers of diameter a 
and separation d, all moving in a fluid with a velocity vo- The 
velocity is measured at a point a distance r from the axis. 



II. ENHANCEMENT OF STREAMING DUE TO 
HYDRODYNAMICS 

Consider impellers to be objects each with maximum 
linear dimension a, and with a mean spacing of d ar- 



ranged in a straight line as shown in Fig [T] These im- 
pellers are pictured as spheres, but hydrodynamics is not 
sensitive to the exact shape of an impeller, as it depends 
mainly on its maximum linear dimension [loj . We will 
now analyze the amount of streaming due to motion of 
these impellers moving on a single microtubule. We ini- 
tially consider the impellers to be much larger than the 
kinesin molecules so that they dominate the hydrody- 
namical response of the fluid. 

If spherical impellers were close-packed along the mi- 
crotubule, that is a w d, then this problem would be 
equivalent to a single rod of length L moving at constant 
velocity Vq in the fluid in a direction parallel to its long 
axis. In that case, the velocity field for distances r L 
has only a weak logarithmic dependence of r, meaning 
that up to a correction of order ln(L/a), the velocity 
field is only weakly dependent on distance and of order 
Vo- (Here we take the velocity of the fluid to go to zero 
far from the rod.) This is related to the well known result 
that the drag on a rod of length L is of the same order as 
that of a sphere of diameter L despite the latter 's much 
greater volume Such a system of densely packed 

impellers would be very efficient at driving fluid motion 
and only a low density of microtubules would be needed 
for cytoplasmic streaming. 

Now consider the more realistic case in which the im- 
peller size is less than the spacing between them, which 
is much less than the length of a microtubule, that is 
a < d <C At a distance r ^ d, the velocity field will 
behave just as in the closed-packed case except appear to 
have a diminished impeller velocity. That is for d <^ r <C 
i, the magnitude of the velocity w(r) = vof{a/d)g{r), 
where g{r) contains all the distance dependence of the 
velocity field, and f{a/d) is how the velocity scales with 
the ratio of a/d. In the limit where a/d is very small, the 
system becomes dilute and the velocity between impellers 
will decay to zero. We recover the motion of isolated im- 
pellers in this case, where it is well known (e.g. Stokes' 
drag) that the velocity field is proportional to a. There- 
fore for small argument x, f{x) is linear (in other words, 
the fluid velocity is proportional to a.) 

Using this general argument we conclude that, inde- 
pendent of the exact shape of the impellers, the flow ve- 
locity is reduced from the closed-packed case by a factor 
^ a/d. The effect of the impellers only starts decreas- 
ing substantially at a distance of order the length of the 
microtubule L, below which it should only have a weak 
logarithmic dependence. In other words, for a spherical 
region just enveloping a microtubule, the fluid velocity is 
slowly varying and reduced from the kinesin motor ve- 
locity by a factor of order a/d. 

Of course there are many microtubules in these cells. 
To understand how this affects the above analysis, con- 
sider all space filled with an infinite forest of them all 
oriented in the same direction. First if we ignore the 
microtubules and just consider the spherical impellers, 
then if we move to a reference frame moving with the 
impeller velocity, the system is static and the velocity 
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everywhere is zero. Therefore in the original reference 
frame, the fluid is also moving uniformly at the impeller 
velocity. This is not correct because we have ignored the 
hydrodynamic drag of the microtubules represented by 
the line going through the spheres in Fig. [1] To esti- 
mate their effect on the fluid velocity, denote the drag 
coefficient on an impeller by 6/ and that of a section of 
microtubule length d by bm- Then we go to a reference 
frame velocity v such that the total force acting on the 
combined system of impellers and microtubules is zero, 
so that bj{vo — v) — Bmv = 0, or v ~ bjVQ/{hi + 6m)- 
Because the net force acting on this system is zero in 
this frame, v is the velocity of the fluid far from the mi- 
crotubules. Because within logarithmic corrections, the 
drag coefficients arc proportional to the maximal linear 
dimensions, then a conservative estimate of this speed is 
of order voa/d. The exact formula depends on the shape 
of the impellers. This argument assumes an infinite vol- 
ume of microtubules but the corrections to this due to 
the finite nature of the system are not important for the 
estimates we are making. 



III. TRAVELLING WAVE INSTABILITY OF 
MICROTUBULES 

We now show how the kinesin generated tangential 
forces on microtubules give rise to travelling wave con- 
formations and calculate their angular and spatial fre- 
quency. A microtubule has a configuration r(s) param- 
eterized by an arclength s and, at long enough length 
scales, can be modeled as being inextensible, that is 
|9r/(?s| = 1 with an elastic bending constant C. The 
inextcnsibility is enforced by a position dependent ten- 
sion T{s). There is also a force acting on the microtubule 
as a result of kinesin walking along it. The magnitude of 
the force is proportional to the local speed and the size of 
the kinesin-driven impeller and the direction of the force 
is tangent to the microtubule, which, we will see, has the 
effect of making it buckle. We denote this with a force 
per unit length of fk- We also include a force due to the 
cytoplasm streaming at a velocity which we take to 
be in the k direction away from the minus end and this 
force tends to straighten the microtubule. This leads to 
the equation 



The identities of impellers in this system are still un- 
known but there are many possible candidates. Anything 
with large linear dimensions in at least one direction will 
give a large hydrodynamic radius a (Toj . The other re- 
quirement is that it can attach to kinesin. In fact it is 
possible that the impellers in this situation could them- 
selves be microtubules that arc not attached to the cor- 
tex [ij, [l^ . Experimental estimates of the cytoplasmic 
streaming velocity are approximately 0.25/im/s, about i 
to ^ the typical velocity of a kinesin molecule. This sug- 
gests that a/d is j or greater. For example, if we take 
the maximum dimensions of an impeller to be 250nTO, 
this predicts a spacing between impellers of Ifim or less. 



Another important biological issue is the necessity to 
have some microtubules in direct physical contact with 
the cortex of the oocyte, for example by tethering or 
by frictional forces. A free floating microtubule with ki- 
nesin moving on it will apply zero net force to the fluid. 
This is a simple consequence of Newton's third law, or 
equivalently, conservation of momentum. This does not 
contradict the fact that bacteria are able to swim: the 
force propelling the bacterium forward is countered by an 
equal and opposite force on the environment, leading to 
velocity fields that are dipolar at large distances. Unlike 
the case analyzed above, this will not lead to long range 
hydrodynamic motion of the fluid and will not lead to 
efficient cytoplasmic streaming by relatively few motor 
proteins. Contact with the cortex is crucial as it allows 
for transfer of force from outside of the oocyte to the cyto- 
plasm enabling fast streaming of the bulk to be powered 
by a much smaller volume of kinesin driven impellers. 
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where j/ is a hydrodynamic drag coefficient per unit 
length. 




FIG. 2. The motion of a microtubule with a constant force 
being applied tangent to its axis. A constant velocity in 
the horizontal direction from left to right has also been ap- 
plied representing the background from cytoplasmic stream- 
ing. The white ball on the left represents the point at which 
the microtubule's minus end is tethered. See supplemental 
movie 1 [l^. In a real oocyte, the radius of curvature would 
correspond to approximately 40/im and the dimensions of the 
entire figure are roughly comparable to the oocyte. The thick- 
ness of the microtubule has been magnified by a factor of 
about 2000 to facilitate viewing. 



We first implemented this equation numerically for a 
range of parameters and enforced boundary conditions 
that tethered the minus end against the cortex, while 
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the plus end was free. Starting from random initial con- 
ditions, the equation rapidly goes to a steady state that 
typically is described by a curve that asymptotically be- 
comes helical for large s and rotates uniformly at con- 
stant angular velocity. The results of a steady state 
configuration are shown in Fig. [2] The supplemental 
movies [l^ discussed below, show microtubule solutions 
for different parameters. Supplemental movie 1 shows 
the full time dependence [l^. The chirality of the helix 
depends on its initial conditions. Therefore the direc- 
tion of rotation is random but stable once steady state is 
reached. 

As the microtubule is made longer, the angular velocity 
and radius of the helix go to a constant limit. However 
the tethered part is not helical but nevertheless rotates 
in synchrony with the rest of the microtubule. These 
dynamics we analyzed in detail (see the supplemental in- 
formation) to find the form of the solution to this equa- 
tion. In particular we show that this equation supports 
travelling waves and find the relationship between the 
angular velocity of rotation uj and the asymptotic radius 
of the helix, and the external velocity fields Vg. In the 
case where Ws = 0, the relationship simplifies to 

uRcu/fk = 1. (2) 

independent of chain length. 

It is interesting to note that for travelling waves, solu- 
tions can be at any scale. There is a continuous family 
of solutions all with the same shape but different scale 
factors. In the case of a helix, many different radii R are 
solutions to these equations. 

What determines the value of R that is selected? As 
with other problems in pattern formation such as the "ge- 
ometric model" or the full dendrite problem [l^, it 
is the boundary conditions that are responsible for the 
unique value of R that is selected. In this case, the 
microtubule minus end is tethered, u(0) = but the 
plus end is free. The solution will only exist for dis- 
crete values oi (i = C / [R? ft)- Numerical analysis gives, 
13 = 0.05 ± 0.0005. This implies that 

R = {C/{Ph)f". (3) 

This analysis was extended to non-zero Vs and as 
shown in the supplemental materials [Toj . does not ap- 
preciably change the estimate we will give below for mi- 
crotubule wave parameters in fast streaming oocytes. 

When the microtubule is tethered to an impenetrable 
surface (in this case, the oocyte cortex) and the external 
cytoplasmic streaming velocity Vs is parallel to that sur- 
face, the form of the solution changes considerably. Nu- 
merical results show that the microtubule becomes com- 
pletely two dimensional, lying close to the surface. For 
low enough Vs , travelling wave solutions are close to pro- 
late cycloids, meaning that the microtubule periodically 
loops back on itself (see supplemental movie 2 [l^ . This 
can be understood analytically. For sufficiently large Vs 
it transitions to other states finally becoming two dimen- 
sional and looking close to a sinusoid (sec supplemental 



movie 3 [16| . This is analyzed in detail in the supplemen- 
tal materials [Tgj . 

We have not included hydrodynamic and stcric inter- 
actions between different microtubules except in the ap- 
proximate way of giving rise to a constant cytoplasmic 
streaming velocity. Given the microtubule density in the 
oocyte, we expect these interactions to be substantial. 
However the wave like motion that we find is quite ro- 
bust. Attaching microtubules together or considering ad- 
ditional forces still leads to periodic or sometimes chaotic 
motion, still at the same characteristic time and spatial 
scales. Therefore despite the simplicity of the model, we 
expect that the basic length and time scales that we pre- 
dict should be quite robust. 



IV. COMPARISON WITH EXPERIMENT AND 
DISCUSSION 

We now check to see whether the above model is con- 
sistent with experimental results. 

Estimates of the microtubule elastic constant C vary 
considerably [l^, [2l| but range mostly within 2 to 4 x 
10~^^A^m^. We estimate the force due to the kincsin per 
unit length, fk to be its velocity times the cytoplas- 
mic viscosity. We will assume as suggested by our above 
analysis, that a/d\s between 1/4 and 1. We will take the 
kinesin velocity Vk to be approximately Ijim/s [TTI. [l^. 
The effective viscosity of the cytoplasm for small par- 
ticles has been studied extensively and appears to vary 
depending on the type of cell on the length scale [1] . Par- 
ticles of different sizes diffusing in the cytoplasm diffuse 
as if the medium had a different viscosity. Its viscoelas- 
tic properties will depend on many factors such as the 
state of gelation of actin filaments [l^] ■ We will assume 
that during fast cytoplasmic streaming, the cytoplasmic 
actin is mainly in a "sol" state allowing fast streaming 
to more readily take place. If this is not the case, the 
effective viscosity could be much higher Q but the dissi- 
pation would increase proportionally requiring a higher 
energy input. For small particles of radius approximately 
20nm, the effective viscosity as measured by diffusion is 
approximately 8 times that of water Q and we will use 
this value for comparison with experiment. This gives 
/j. {a/d)vr] in the range 2 to 8 x 10~^A^/m. Plug- 
ging these numbers into Eq. [3] gives R in the range 30 
to 50/im. We can estimate the angular velocity using 
Eq. [2] (which ignores cytoplasmic streaming). To get 
the largest range of times, we assume that the hydrody- 
namic drag coefficient is due solely to the impellers, so 
that v = {a/d)T]. This is equivalent to w = Vk/R or a 
period of T = 2iTR/vk which ranges from 190 to 310s. 

Serbus and colleagues have observed microtubule be- 
havior in living Drosophila oocytes using GFP-tubulin 
and confocal fiuorescence microscopy [l[. Time-lapse 
movies showed bright fibrous fiuorescence, representing 
groups of microtubules, in a background of diffuse fluo- 
rescence from non-polymerized GFP-tubulin. The micro- 
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tubule patterns changed over time, as did the patterns 
of motion of organelles that either excluded the GFP- 
tubulin or were themselves auto- fluorescent. Bends in 
the microtubules were at times visible in the x — y opti- 
cal plane and some remained visible within that plane for 
30 — 90s, allowing measurement of a radius of curvature 
that wc expect to approximately correspond to the R in 
the above analysis. R was measured to be 19.5 ± 6.5/im 
with 9 measurements (standard error 2.1fim), and the 
wave velocity was v = 0.265 ± 0.04/im/s with 4 measure- 
ments. Because these waves appear roughly sinusoidal, 
we can also estimate the period T. Assuming a wave- 
length of 4i? then the measured velocity gives a period of 
T = 294s with an even larger error bar considering the 
fact that this assumes a waveform that is probably not 
accurate. Nevertheless it suggests a characteristic time. 

The motion studied here can be contrasted with cil- 
iary motion (23j which has a period on the order of .05s. 
Clearly the two mechanisms have fundamentally differ- 
ent explanations. The experimental time scale points 
to the mechanism described here rather than ciliary mo- 
tion. The agreement found with experiment in our above 
analysis is to some extent fortuitous, given the large un- 
certainties in the experimental system. But nevertheless, 
it provides evidence that the simple mechanism proposed 
is the origin of the behavior seen in cytoplasmic stream- 
ing. The experimental evidence [l| that dynein inhibits 
streaming and that minus ends of microtubules are in 
contact with the cortex [l^, both support this hypothe- 
sis as well. 

More generally, the above analysis elucidates, at a 
qualitative level, the phenomena seen in experiment [J, 
of strong cytoplasmic streaming occurring concomitantly 
with wave-like motion of microtubules. The movement 
of a relatively few number of kinesin motors on a micro- 
tubule couples to the fluid causing a bulk hydrodynamic 
flow. The force a kinesin exerts on a microtubule sets the 
microtubule into motion causing it to execute wave-like 
motion. This in turn causes the flow lines around micro- 
tubules to be time dependent, making the microtubules 
act as stirrers for the surrounding fluid, which leads to 
chaotic flows and strong mixing 0, Q . 




FIG. 3. The configuration a rotating microtubule in a gliding 
assay (2^, [2^ with the leading end stuck to a kinesin protein 
is shown in (a), (b) Shows the two dimensional simulation of 
Eq. [T] showing a similar shape. 



Can such motion be directly observed in an in vitro 
experiment? In fact it has been seen frequently as an un- 



wanted artifact of sample preparation. Microtubule glid- 
ing assays were pioneered more than two decades ago [l^l 
and have become a standard technique to better under- 
stand the motion of kinesin motors. When a solution mi- 
crotubules and ATP is placed on a glass surface on which 
there is a high density of kinesin molecules, the kinesins 
propel the microtubules with their minus ends leading. 
Occasionally the minus end sticks, perhaps to damaged 
kinesins, and continued forces of the trailing portion of 
the microtubule cause flexing and curve generation. 

Microtubule gliding assay videos of high quality 
provide an excellent source of data for investigating the 
two dimensional version of the problem studied here [26| . 
The microtubule is being pushed by kinesin molecules 
that on average, apply a force tangent to the microtubule. 
When an end is tethered by sticking to a kinesin molecule, 
as noted above, this leads to a physical situation that is 
well described by Eq. [TJ It should be emphasized that 
the values of the tangential force fk and the drag co- 
efficient 1/ will be much larger than in the oocyte case. 
Consequently, sometimes microtubules flex,, forming a 
G-shaped curve that rotates around the anchored mi- 
nus end of the microtubule. Fig. |3] shows a comparison 
of snapshots of the microtubule with that of a simula- 
tion of Eq. [T] The resemblance is quite striking and 
provides evidence that the theoretical modelling of this 
bending instability is valid. Note that the shape obtained 
is independent of the values of model parameters. From 
measuring the radius R of the G-like configuration and 
its angular velocity of rotation, we can compute both fk 
and I' from Eqs. 13] and [5J giving v = ujC/{f3R'^) and 
fk = lyRuJ. Measurements give R w 1.33 ± 0.2/iTO and 
the period T = 12 ± Is. Therefore ly = 133Ns/m'^ and 
/fc = 9.3 X 10^^ N/m, both estimates are accurate within 
a factor of 2 assuming that C has been well determined. 
The force exerted by a kinesin motor in this situation is 
close to the experimentally determined stall value of ap- 
proximately fm = 5pN [IJl. Assuming that each motor 
attached can exert this force, this gives an average dis- 
tance between motors of fm/ fk = 53nm. This assumes 
that the attachment of the motors is flexible enough to 
allow microtubule interactions capable of exerting this 
stall force. If this is not the case, we expect the kinesins 
to be spaced more closely but by no more than half this 
amount. This distance seems reasonable given the high 
density of kinesin molecules attached to the glass, and the 
size of the a — (3 tubulin dimer which is approximately 
8nm in length. Thus the above considerations can give 
a quantitative measure of the forces and drag in these 
gliding assays. 

Nevertheless, it would also be of interest to verify the 
travelling wave solutions experimentally in three dimen- 
sions for isolated microtubules and kinesin, for example, 
by utilizing a single molecule optical trap, or by some 
other means, allowing for a more rigorous test of the pre- 
dictions made here. 

It is of interest to speculate on the reason for this mi- 
croscopic mixing mechanism. If the microtubules were 
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attached from the plus end, there would be no dynamical 
instability from kinesin generated forces but the cyto- 
plasmic streaming would be more efficient, with greater 
transfer of force from outside the oocyte, because the 
tension at the tether point would increase linearly with 
the length of the microtubule, whereas the tension sat- 
urates to a constant value in the case considered above. 
However as discussed earlier, complex stirring patterns 
are expected to be far more efficacious in mixing than 
steady-state flows. This agitation could then serve an 
important biological function. 

In summary, we have found a very simple set of con- 
ditions for generating efficient mixing in large eukaryotic 
cells. If microtubules are aligned so that their minus ends 
are in contact with the cell perimeter and dynein is sup- 
pressed, this will lead (a) to cytoplasmic streaming in 
the bulk of the cell and (b) wave-like motion of micro- 
tubules. The latter will promote efficient mixing of the 
cell's contents by inducing chaotic flows. The efficiency 
of the cytoplasmic streaming is due to the long range 
nature of hydrodynamic coupling and the indirect link- 
age of microtubules to the oocyte surroundings via the 
cortex, leading to momentum transfer. The wave-like 
motion is due to an instability in the dynamics which 
breaks chiral symmetry and is a consequence of the tan- 
gential forces exerted by motor proteins on the elastic 
microtubules. Given the uncertainties in the experimen- 
tal system and approximations made in the analysis. The 



wave- like motion seen in experiments agrees well with our 
theoretical predictions suggesting these effects represent 
a robust phenomenon. We are now engaged in experi- 
mentation to gather more accurate data on microtubule 
behavior in vivo and in vitro under various conditions 
in order to further understand the mechanisms for fast 
cytoplasmic streaming discussed here. 
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Note: After the completion of this work, we came 
across work of Bourdieu et al. [2^ which analyzed their 
experimental results on spirals in mysosin and kinesin 
gliding assays using the two dimensional version of Eq. 
[TJ obtaining the identical scaling, and simulating it us- 
ing a nearly identical approach. This provides further 
evidence that the mechanism we have proposed for cyto- 
plasmic streaming is physically viable. 
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